BUTP-95/4 



Chiral symmetry and spectrum of Euclidean Dirac 

operator in QCD Q 



A.V. Smilga 

Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, Bern 

CH-3012, Switzerland^ 

March 1995 



Abstract 

Some exact relations for the spectral density p{X) of the Euclidean Dirac 
operator in QCD are derived. They follow directly from the chiral symmetry 
of the QCD lagrangian with massless quarks. New results are obtained both 
in thermodynamic limit when the Euclidean volume V is sent to infinity and 
also in the theory defined in finite volume where the spectrum is discrete and a 
nontrivial information on p(A) in the region A ~ 1/(| < qq >o \V) (the charac- 
teristic level spacing) can be obtained. These exact results should be confronted 
with "experimental" numerical simulations on the lattices and in some partic- 
ular models for QCD vacuum structure and may serve as a nontrivial test of 
the validity of these simulations. 



1 Introduction 

The notion of spectral density of the Euclidean Dirac operator in QCD has been 
brought into discussion some years ago in the pioneering paper by Banks and Casher 
]1|]. They have got the famous formula 

p(0) = --<qq > (1.1) 

7T 
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relating the infrared limit of the averaged spectral density with the quark condensate 
- the order parameter for the spontaneous chiral symmetry breaking. p(A) is defined 
as 

p(\) = ^<uj(\,A)> A (1.2) 

where 

u(X,A)= y £S(X-X n ) (1-3) 

n 

is the microscopic spectral density of the Dirac operator in a given background field 
A in a finite Euclidean volume V. The averaging is taken over all gluon fields with 
the standard Yang-Mills measure involving the determinant factor for Nf fermions 
with the common mass m. The combination Vp(X)dX defines the average number of 
eigenvalues of Dirac operator in the interval (A, A + dX). 

Let us recall how the Banks-Casher relation fll.lj ) (which holds in the thermody- 
namic limit V — > oo, with the mass m being kept small but fixed) is derived. Treating 
the gauge field AHx) as an external field, the fermion Green's function is given by 

S A (x, y) =< q(x)q(y) > A = £ (L4 ) 

V m - iX n 

where u n (x) and X n are eigenf unctions and eigenvalues of the Euclidean Dirac oper- 
ator: 

Tpu n {x) = X n u n (x) (1.5) 
Note that the spectrum of the theory enjoys the chiral symmetry: 

Ipn -> 7 5 ^n, A n -> -A n (1.6) 

And, except for zero modes, the eigenfunctions occur in pairs with opposite eigenval- 



ues. Setting x = y and integrating over x, the representation ( |1.4|) therefore implies 
If 2m „ 1 

— / dx < q{x)q(y) > A = -— ^ - 1.7 
V Jv V C^ m 2 + X 2 n 

where the zero mode contributions have been dropped out (it is justified in the ther- 
modynamic limit. See Ref.0 for details). In the limit V — > oo, the level spectrum 
becomes dense and we can trade the sum in the r.h.s. of Eq. (|1.7|) for the integral. We 

get 

< qq > = -2m f°° dX ^ (1.8) 
Jo m z + X z 

The integral ( pL8|) diverges at large A. For A ^> Aq C d the spectral density is 
not sensitive to gluon vacuum fluctuations and behaves in the same way as for free 
fermions: 

p free {X) „ ^X_ ^ 
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This perturbative quadratically divergent piece in < qq > is proportional to the 
quark mass and is related just to the fact that mass terms in the lagrangian break the 
chiral symmetry explicitly. To get the truly non-perturbative quark condensate which 
is the order parameter of the spontaneous chiral symmetry breaking, this perturbative 
divergent part should be subtracted. As a result, the quark condensate is related to 
the region of small A: A ~ m <C ^qcd- It is not difficult to see that the finite 
mass-independent contribution appears if p(0) 7^ and the relation holds. 



2 Thermodynamic limit: A ^ 0. 

In this section, we derive a new formula describing the behavior of the spectral density 
p(A) in the region m C A < ^qcd in the thermodynamic limit V — > oo ||. To this 
end, let us consider the 2-point correlator 

K ab = Jd*xJ d 4 y < 0\S a {x)S b {y)\0 > (2.1) 

with S a = qt a q where t a is the generator of the flavour SU{Nf) group. The correlator 
K ab can be evaluated in the same way as < qq >o- First, fix a particular gluon 
background. As Tr{t a } = 0, only the connected part depicted in Fig. la contributes, 
and one obtains 

Kf = -J d 4 x J d*yTr{t a g A (x,y)t b g A (y,x)} (2.2) 

Substituting here the spectral decomposition for the Green's function ( |1.4[ ) and taking 
into account the chiral symmetry of the spectrum (|1.6| ), we get 



K a A b = S ab E /"I,' , A 1 (2-3) 



Rewriting the sum as an integral as in Eq.( |1.8| ) and averaging over the gluon back- 
ground, one gets 

1 rsab X ab rp(X)(m 2 -X 2 ' 



On the other hand, the correlator ( |2.1| ) of the colorless scalar currents can be evaluated 
by inserting the complete set of physical states. For small m, i.e. close to the chiral 
limit, a distinguished position among the latter belongs to the Goldstone states which 
appear due to spontaneous breaking of the chiral symmetry of the QCD lagrangian. 
Low energy properties of Goldstones are fixed by chiral symmetry, and some exact 
calculations are possible. The chiral lagrangian has the form Q 
1. 



Ti{(d^U ] )(d^U)} + B Re Ti{MU^} 



+ higher order terms (2.5) 
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where B = — < qq >q F~ 2 , U = exp{2i<p a t a / F n } , and Ai is the quark mass matrix. 
If M. = diag(m, . . . , m), the lagrangian ( |2.5|) describes Nj — 1 (quasi-) Goldstone 
states with the common mass 

Mj = 2mB (2.6) 

Consider now the graph in Fig. lb with 2-goldstone intermediate state contributing 
to the correlator (|2.4j) (obviously, one-goldstone state does not contribute since 
< 0\S a \<p b >= 0). To calculate it, we need to know the vertex < 0\S a \<ft b <j) c >. It can 
be easily found from the generating functional involving external scalar sources. The 
latter is obtained substituting in the effective lagrangian Q2.5|) the mass matrix M. by 
A4 + u a t a where u a is the source for the scalar current S a . In this way, one gets 

< 0\S a \<f) b <f) c >= Bd abc (2.7) 

It is very important that the vertex is nonzero only for three or more flavors. Now 
we can calculate the graph in Fig. lb absorbing its ultraviolet divergence into local 
counterterms contained in higher-order terms in the effective lagrangian Q2.5|) [[| . The 
result reads 

= - ^^y ** (2-8) 

For massless goldstones, the graph in Fig. lb does exhibit a logarithmic infrared sin- 
gularity reflected in the factor lnM? on the r.h.s. of Eq. ( |2.8|) . The same circumstance 
makes our calculation self-consistent: since infrared singularity arises from the low 
momenta region, higher derivative terms in Eq.( |2.5| ) can be neglected. 

Now, let us compare Eq. (|2.4j ) with Eq.( |2.8|) . Note first of all that, in contrast 
to the integral in Eq.( |l~8|) , the constant part p(0) does not contribute here at all: 



ft^^O (2.9) 

Thus, only the difference p(A) — p(0) is relevant. It is easy to see that, in order to 
reproduce the singularity oc In M| oc In m, one should have 

p(X)-p(0) = CX (2.10) 

at small A. Comparing the coefficients at the logarithms mm and In Mj in Eq. ( |2.4j ) 
and Eq.( [2.8|) , we arrive at the result 

The structure |A| appears due to the relation p{— A) = p(A) being implied by the chiral 
symmetry (|1.6j). Note again that this non-analyticity does not appear at Nf = 2. 
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Recently, the preliminary, not yet published data for the spectral density p(A) 
calculated on the lattice with two dynamical fermion flavours appeared || . The same 
quantity has been also evaluated in the instanton-antiinstanton liquid model for QCD 
vacuum configurations |J . Results of these two numerical studies are completely dif- 
ferent. Lattice calculation failed to reproduce the exact QCD relation (|2.11|) (they 
exhibit a large nonzero slope for p(A) which should be absent for Nf = 2). Probably, 
it is due to the fact that , in the region of lattice parameters used in [||], the thermo- 
dynamic limit was not reached yet and the finite volume effects were still important 
0. Obviously, further studies in this direction are highly desirable. 

On the other hand, the calculations in the instanton model agree well with 
theoretical predictions. The results for the spectral density obtained in || for different 
number of light quark flavors are presented in Fig.2. Before comparing them with the 
theory, two remarks are in order. 

1. Only non-perturbative instanton-driven part of spectral density has been deter- 
mined in ||. So, the perturbative effects which give the dominant contribution 
to the spectral density at large A were not taken into account, and the compar- 
ison with theory makes sense only in the small A region. 

2. The abrubt falloff and vanishing of p(A) at A = as measured in |J is the finite 
volume effect . The comparison should be done in the region of not yet too small 
A. 

From the graphs in Fig.2, one sees that, for Nf = 2, the spectral density is 
practically flat, while for Nf = 3 it rises with a nonzero slope. This is exactly what 
the formula fl2.11|) requires. 



The derivation of the result ( [2.11D assumed the spontaneous breaking of chiral 



symmetry and the existence of the Goldstone bosons. Thus, it refers only to the case 
Nf > 2. For Nf = 1, there is no theoretical result, and no comparison can be done. 
It is amusing, however, that, if trying to "continue analytically" the formula 



down to the point Nf = 1, one would get a negative slope for p(A) which agrees again 
with the instanton simulations. 

Qualitatively, the different A - dependence for different Nf is natural. Mean 
spectral density p(A) is obtained after averaging of the microscopic spectral density 
( |1.3|) over gluon fields with the weight which involves also the quark determinant 



[det(z^)] 



Nt 

N f 



n \ a 

A„>0 



(2.12) 
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The larger Nf is, the more the region of small A is suppressed. There are good reasons 
to expect that, at Nf = 0, with no suppression at all, quark condensate = p(0) is 
infinite in the thermodynamic limit M. Q 



3 Finite volume: the partition function. 

In the rest of this talk, I derive and discuss some relations for the spectral density in 
the region of very small A ~ 1/| < qq > \V. These relations refer not to physical QCD 
in the infinite volume, but to the theory defined in the finite Euclidean box. But, as 
our main goal is to discover islands of firm theoretical ground in the foreboding sea 
of numerical simulations , and the latter are done exclusively in finite volume, these 
intrinsically finite volume results can serve this purpose in exactly the same way as 
the result ( |2.11| ) derived in the infinite volume. 

But before proceeding to the Dirac operator spectrum, we are in a position to 
study the more basic quantity — namely, the QCD partition function at finite volume. 
Specifically, we will be interested with the dependence of the partition function on 
the quark mass matrix Ai and the vacuum angle 9. 

Consider first the theory with only one light quark flavor. It involves a gap 
in the spectrum {U(l) axial symmetry is broken not spontaneously but explicitly by 
anomaly, and no goldstones appear). Thus, the extensive property for the partition 
function holds: 

Z ~ exp{-e„ ac (m, 9)V} (3.1) 

If L Aq^ d , the finite volume effects in e vac are exponentially small ||. Ward 
identities dictate that e vac can depend on m and 9 not in an arbitrary way, but only 
as a function of a particular combination me 10 . Expanding e vac (me l6 ) in Taylor series 
(that makes sense as long as m < Aqcd) and bearing in mind that e vac is real, we 
get 

1 r^j gSVmcose + 0(m 2 ) ^ 2) 

The parameter £ is nothing else but the quark condensate (up to a sign). Really, 

< q q > = - J- In Z = -S (3.3) 
dm 

Consider now the same problem but for several quark flavors. Then spontaneous 
breaking of chiral symmetry occurs, goldstones appear, there is no gap in the spectrum 



3 The instanton calculations of Rcf.|| do not exhibit the infinite quark condensate for the 
quenched theory, but only a much more dramatic rise of the spectral density in the low A region. 
But, probably, the finite value for p(0) is, again, a finite-volume effect. 
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, and the extensive property ( |3.1| ) does not hold anymore. More exactly, if quark 
masses are non-zero, the property (3T) still holds when the length of the box L is much 



larger than the Compton wavelength of Goldstone particles ~ ~ 1/ ymAgcD- 

We will be interested, however, with the intermediate region 

Aqcd < L < - ] (3.4) 
^mA Q cD 

where the effects due to goldstones on the volume and mass dependence of the par- 
tition function are crucial. 

Fortunately, the goldstone properties at low energies (and, in the region L ^> 
Aqc D , only the low-energy properties of the Goldstone fields are relevant) are well 
known. They are described by the effective lagrangian (|2.5| ) where, at nonzero 9 , M. 
should be substituted by Me i9/N f ( again, the occurence of this particular combination 
is dictated by Ward identities). The partition function is given by the functional 
integral 

[dU] exp j-i^ 2 J d 4 x -Tr-T^C^^C/)} + B Re Tr{Me i0/N ftf j (3.5) 
The crucial observation is that, as long as L <C — ~ , = , only zero Fourier 

6 m <t> y/mAgco ' J 

harmonics Uq of the field U(x) is relevant, and the contribution of the higher har- 
monics in the integral ( |3.5| ) is suppressed. Then the functional integral is transformed 
into the finite-dimensional integral 

Z~[ dfj,{U Q ) exp{VE Re Tt {Me ie/N fU^}} (3.6) 

JSU(N f ) 

The integral extends over the group SU(Nf), and dfi(Uo) is the corresponding Haar 
measure. For Nj = 2, the integral ( |3.6|) can be done explicitly. For Nf > 3, it is not 
possible to do analytically in the general case. Sometimes (for the case M. = ml), 
certain beautiful results can be obtained 0, but, to derive the sum rules for the 
small eigenvalues of the Dirac operator we are after in this talk, the closed analytic 
expression for Z(V,A4,9) is not actually required. 

The partition function ( |3.6| ) is a periodic function of 9 with the period 2ir. Thus, 
it can be presented ClS db Fourier series 

oo 

Z{9)= £ Z v e lvB (3.7) 

v=— oo 

so that 

Z u = — d9e" iv6 Z{9)d9 ~ / ddU) (det UY exp{VZ Re Tr{MUH} (3.8) 

2,71 JO JU(Nf) 



with U = Uoe- ie ' N f. 
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4 Finite volume: sum rules. 



The main idea is basically the same as in Sect. 2 — to compare the expression (|3.8| ) 
derived in the chiral theory with the quark-gluon representation of the same quantity. 
Z v may be writen as a functional integral over the quark and gluon fields restricted 
on the topological class with a given winding number 

" = ^2 / G luG% (4.1) 

The integral over the Fermi fields produces the determinant of the Dirac operator, 
and we get for v > 

Z u = [ [dA}e-H d4xG ^ G ^(det f M) u JJ det f (\ 2 n + MM ] ) (4.2) 

A n >0 

where the factor (detfAi) u arises due to the fermion zero modes X n = which appear 
on a topologically non-trivial background due to the index theorem. (If v is negative, 
the factor [detfM) v is to be replaced by (detfj\4^)~ u ). 

Let us expand now Eqs.Q3.8p and ( |4.2j ) in quark mass and compare the co- 
efficients of (detfA4) v Tr{M.A4'}. On the chiral side of the equality, we get some 
group integral which can be done explicitly. On the quark and gluon side, we get the 
expression involving < J2x n >o V^n > f where the average is defined as 

I[dA]fe~ l J^o^ {Uxn>Q Xlff 
V f[dA)e-U^c^ {Il , n> oK) Nf 

where the path integral is done over all gauge fields with the topological charge v. 
Skipping the technical details, we present the final result for the simplest sum rule 
thus derived: 

(VIL^ (4.4) 

where k = | v \ + Nf . 

Expanding the expressions ( p.8| ) and Q4.2| ) further in quark mass and comparing 
the coefficients of the two group invariants (det f M) u (Ti{M M ] }f and 
{det f M) u Tr{MM ] MM ] } , we get two different sum rules for the inverse fourth 
powers of eigenvalues 

(a5o k) v = mk 2 -iy ( ( A 5o ^) )„ = i6(* 2 -i) (4,5) 

(Two different sum rules ( |4.5[ ) are obtained when Nf > 2. If Nf = 1, there is 
no difference between (Tr{.M.A/^}) 2 and Tr{.M.M^ M.M^} and only one sum rule 
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involving inverse fourth powers of A n for the difference of the two sums entering 



Eq.(4.5) can be derived). Further expansion in A4 provides the sum rules with inverse 
sixth powers of eigenvalues etc. 

Note that the sums Q4.4| ) , (|4.5| ) are saturated by few first eigenvalues. There- 
fore the sum rules provide the information on the very bottom of the finite-volume 
spectrum. 

The theoretical results ( |4.4|) , ( |4.5| ) have not yet been confronted with lattice 
experiments. However, the numerical calculations in the instanton model in the case 
v = with different number of quark flavors are available. They are in a good 
agreement with the theory (see Fig.3 and the discussion thereafter). 



5 Stochastic matrices 

The sum rules ( |4.4|) , (|4.5|) are exact theoretical results. They follow just from the form 
of QCD lagrangian and the assumption (which is the experimental fact in QCD) that 
the spontaneous breaking of chiral symmetry occurs. It is interesting, however, that 
much reacher information on the spectrum of Dirac operator in the region A ~ 1/ {^V) 
can be obtain if invoking a certain additional assumption on the properties of the 
functional integral in QCD ||10|| . 

Consider the Dirac operator in a particular gauge field background with the 
topological charge v. For simplicity, let us assume that the mass matrix is diagonal 
and real. Let us choose a particular basis £ = {ip n {x)} and present the Dirac operator 
— % Tp+mf as a matrix in this basis. The particular choice of the basis is not important 
for us, but a reader familiar with the model of instanton-antiinstanton liquid may 
think of ip n {x) as of zero-mode solutions for individual instantons and antiinstantons. 
What is important is that the basis S is chosen in such a way that all states ip n (x) 
have a definite (left or right) chirality. In that case, the number of left-handed states 
rii and the number of right-handed states tir in the basis should be related: 

n L - n R = v (5.1) 

Thereby, the existence of exactly v left-handed [y > 0) or —u right-handed [y < 
) zero eigenvalues of the massless Dirac operator —% lp is assured. In this basis, the 
full Dirac operator —iT/> + nif can be presented as a matrix 

M I nit iT \ , . 

-iV + m, = f (5.2) 

V zTT m f J 
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where T is a general rectangular complex matrix x ur. We assume ul^Ur to be 
very large and will eventually send them to infinity. In this limit, the basis £ spans 
the whole Hilbert space and the representation ( |5.2| ) is just exact. 

The partition function involves the Dirac determinant and the integral over all 
gluon configurations. If the basis £ is fixed (the same for all gluon field configurations), 
the path integral over the latter can be traded for an integral over the matrices T. 
Thus, we have 

Z v = [vTP(T)Jldet( mf %T ) (5.3) 
J t \ iTi m f J 

, the condition ( |5.1| ) for the dimensions of the matrix T being assumed. 

And now comes the crucial assumption. Let us assume that the measure P(T) 
has the simple form 

P{T) = exp |-^Tr{TT f } + o(TT^ (5.4) 

The constant A will at the end determine the spectral density p(0) in the thermody- 
namic limit. 

The form ( |5.4|) of the weight is rather natural. It respects chiral symmetry, 
analytic in T , and is conceptially similar to the Gibbs statistical distribution. How- 
ever, a rigourous proof that P(T) should have this form is absent by now. Therefore, 
the results for the microscopic spectral density p(A ~ \/Y>V) derived in (|[0|) which 
follow from the representation ( |5.3j ) for the partition function and the heuristic as- 
sumption ( |5.4j ) have not the same status as the sum rules ( |4.4|) , ( |4.5| ) which are the 
theorems of QCD. 

Let us explain (very sketchy) how the results for the microscopic spectral density 
have been derived. For simplicity, assume that v = so that T is a square matrix: 
n L — n R = n - P{T) and the Dirac determinant depend only on the eigenvalues A n of 
the matrix T. Therefore, it is convenient to write the integral over T>T as the integral 
over eigenvalues and do the integral over remaining angular variables on which the 
integrand does not depend. We have Q 

VT = C - A; 2 ) 2 J] hd\ k (5.5) 

k<l k 

4 This transformation and the whole stochastic matrix technique is well known for nuclear physi- 
cists and solid state physicists who are concerned with the problem of quantum dynamic for stochas- 
tic hamiltonia. The particular problem studied here coincides with the known problem of " Gaussian 
unitary ensemble" . The vanishing of the Jacobian at Afc = A; reflects the notorious repulsion of the 
levels. For an extensive review of the relevant mathematics sec [01. 
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Thus, joint eigenvalue distribution is just 

A n 

p(A 1; . . . , A n ) ~ U(Xt - A*) 2 II A fc exp{-- £ A 2 } U(\t + mj) (5.6) 

k<l k 71 k=l / 



The microscopic spectral density is obtained after integration of (|5.6|) over all eigen- 
values but one. After some work, one gets the answer in the limit n — > oo. In |10| 



the result was obtained in the particular case v = 0. Being expressed via physical 
variables V (the volume of the system) and £ (the quark condensate), A = (T,V) 2 /A, 
it has the form 

p(A) = — ^— [%r f (pv\) - J Nf+1 (ZV\)J Nf - 1 (EV\)\ (5.7) 

where J^(x) are the Bessel functions. The result (|5.7|) has been obtained in [|1(| under 
the assumption v = 0. In fl2| it has been generalized to the case v ^ 0. [the only 
change is that Nf should be substituted by Nf + \u\ as in Eqs. ( |4.4j , 

At very small A l/VH, the spectral density is suppressed p(X) ~ \ 2N f +1 (the 
suppression is due to the determinant factor Q2.12D which punishes small eigenvalues) . 
But then it rises and after several oscillations with decreasing amplitude levels off at a 
constant p(0) = S/7r as it should in the thermodynamic limit V — > oo. The function 
( |5.7|) with Nf = 0, 1,2 together with the results of numerical instanton calculations 
are plotted in Fig. 3 borrowed from Ref.||. Likewise, one can integrate Eq.( |5.6j ) over 
all eigenvalues but two and derive the expression for the correlation function p(Xi, A2) 
etc. 

The sum rules derived in the previous section can be expressed as the integrals 
of p(A), p(Ai, A2) etc. with a proper weight. For example, the simplest sum rule is 

'Jo -J^-^Nj (5 ' 8) 

which coincides with ( |4.4| ). Thus, the excellent agreement of the model instanton 
calculations with the stochastic matrix model results displayed in Fig. 3 implies also 
the agreement of the former with the exact theoretical results ( [4.4|) , (|4.5|) etc. 




6 Exotic theories. 

The same program as in QCD can be carried out in other gauge theories with non- 
standard fermion content. The results derived earlier depended on the assumption of 
the standard pattern of chiral symmetry breaking 

SU L (N f ) ® SU R (N f ) -> SU v (N f ) (6.1) 
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so that N* — 1 Goldstone bosons living on the coset, which is also SU(Nf), appear. 
However, the breaking according to this scheme only occurs for fermions belonging to 
the complex (e.g. fundamental) representation of the gauge group. The SU(N C ) group 
involves truly complex representations when N c > 3. For N c = 2, the fundamental 
representation (and also other representations with half-integer isospin) is pseudoreal: 
quarks and antiquarks transform in the same way under the action of the gauge group 
and the pattern of chiral symmetry breaking is different leading to different sum rules. 
A third pattern of chiral symmetry breaking is for fermions in the real (e.g. adjoint) 
representation leading to yet another class of sum rules. 

Consider first the case of pseudoreal fermions. It may have a considerable 
practical importance as the numerical calculations with the SU(2) gauge group are 
simpler than with 577(3) group, and it may be easier to confront the sum rules with 
lattice simulations. The true chiral symmetry group of the QCD-like lagrangian but 
with SU{2) gauge group is SU(2N f ) rather than just SU L (N f ) ® SU(N R ) ® U v (l) 
(it involves also the transformations that mix quarks with antiquarks). The pattern 
of spontaneous symmetry breaking due to formation of quark condensate is [pl| 



SU(2N f ) -> Sp(2N f ) (6.2) 

The breaking (|6.2j ) leads to the appearance of 

(2N f ) 2 - 1 - (2N 2 f + N f ) =2N 2 f -N f -l 

Goldstone bosons which are parametrized by the coset SU(2Nf)/ Sp(2Nf). For Nf = 
2 , we have 5 instead of the usual 3 Goldstone bosons. Let us now allow for small 
non-zero quark masses (which break the chiral symmetry explicitly and give non-zero 
masses to Goldstone particles) and study the dependence of the partition function in 
the sector with a given topological charge v on the quark mass matrix M. and the 
volume V. The analog of the result ( |3. 8| ) for the pseudoreal case which is valid in the 
region ( |3.4| ) is the following (see for more details) 



I du(U)(detU)- u exp{ ReTi{MUIU T }\ (6.3) 

JU(2Nf) I 2 J 



where 




(6.4) 



is the symplectic 2Nf ® 2Nf matrix. Formally, the representation ( |6.3|) involves the 
integral over U(2Nf), but actually it is the integral over the coset SU (2N /) / Sp(2N /) 
involving less number of parameters (and the integral over the U(l) part of U(2Nf) 
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which is nothing else but a vacuum angle 9) — multiplying U by a simplectic matrix 
G Sp(2Nf) leaves the integrand invariant. 

Expanding ( |6.3| ) over the group invariants involving the quark matrix M. and 
expanding also the quark-gluon representation of the same partition function, the 
sum rules for the inverse powers of the eigenvalues can be derived. The simplest sum 
rule is 

CSolOtM + sVi) (6 - 5) 

Consider now the third non-trivial class of theories where the fermions belong to 
the real representation of the gauge group. The supersymmetric Yang-Mills theories 
involving Majorana fermions in the adjoint color representation belong to this class. 
We will not assume, however, that the theory is supersymmetric and consider a Yang- 
Mills theory coupled to Nf different Majorana adjoint fermion fields. Also, we will 
assume that the fermion condensate is formed and the chiral symmetry which the 
lagrangian enjoys for Nf > 2 is spontaneously broken Q. 

The chiral symmetry group of such a theory is SU(Nf). Formation of the 
condensate breaks it down to SO(Nf), and 

^2 _ i _ N f (N f - 1) = N f (N f + l) _ 
f 2 2 

Goldstone bosons appear. Again, the partition function Z v is given by the integral 

over the coset SU(Nf)/SO(Nf) which has the form 

Z v ~ [ d/j(U)(detU)- 2vNc exp \VI: Re Tr{MUU T }\ (6.6) 

JU(Nf) 1 > 

The essential difference with the case of fundamental fermions is that the admissible 
values for the topological charge v are not integer but integer multiples of 1/N C . That 
means, in particular, that the function Z(A4,6) given by the Fourier sum ( |3.7] ) is a 
periodic function of 9 with the period 2irN c (not 2tt as before). See Ref.[0 for a 
detailed discussion of this important issue. 

To derive the sum rules, we have to expand ( |6.6j ) over Ai and to compare it with 
the expansion of the quark-gluon path integral for the partition function. And here 
we meet a known problem characteristic for theories with Majorana particles. The 
matter is that the reality condition for the Fermi fields can be fulfilled in Minkowski 



space but not in Euclidean space - it immediately leads to a contradiction [^q] 

5 This dynamic assumption is not so innocent. For example, in N — 2 supersymmetric Yang- 
Mills theory it does not happen [ p^[ . But this theory involves also massless adjoint scalar fields with 
particular interaction vertices. It is quite conceivable that the non-supersymmetric theory with only 
gauge fields and two Majorana adjoint flavors involves chiral symmetry breaking. All the results are 
derived under this assumption. 
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However, the Euclidean path integral for the partition function in a theory with 
Majorana fermions still can be defined by analytic continuation of the Minkowski 
space path integral. In the Minkowski space, the path integral over fermion fields 
is equal to the square root of the corresponding Dirac determinant. The latter can 
be continued to Euclidean space without problems, and the square root is also taken 
easily here due to a notable fact that the spectrum of the Euclidean Dirac operator 
for the adjoint fermions is doubly degenerate: for any eigenfunction u n (x) with an 
eigenvalue A n , the function C~ x u* n (x) (C is the charge conjugation matrix) is linearly 
independent from u n (x) and has exactly the same eigenvalue |§. Thus, square root 
of the Dirac determinant is just 

[det(-ip + M} 1/2 = (det f M) uNc Y[ det f (X 2 n + MM ] ) (6.7) 

A„>0 

where M. = |(1 — 7 5 ).M + |(1 + 7 5 )-Mt, the product counts only one eigenvalue of 
each degenerate pair, and the factor (detfM) uNc reflects the presence of vN c pairs 
of zero eigenmodes of adjoint Dirac operator in the gauge field background with the 
topological charge v |. 

The final form of the simplest sum rule for the eigenvalues of the Dirac adjoint 
determinant is 

~~ 4(\u\N c+ (N f + l)/2) (6 ' 8) 

where, again, only one eigenvalue of each degenerate pair is taken into account in the 
sum. 

The sum rules (|4.4j), (|6.5p, and ( |6.8D can be written universally as 

V 2 Y? 





4{>| + [dim (coset) + 1]/N f } 



(6.9) 



with the rescaling v — > uN c and counting in the sum only one eigenvalue of each 
degenerate pair in the adjoint case. 

The sum rules ( |6.5|) , and ( |6.8|) and their analogs with higher inverse powers 
of A n can be derived also in the stochastic matrix technique, and not only the sum 
rules, but also the expressions for the microscopic spectral density p(A), correlators 
p(A 1; A 2 ) etc. @ 



Mathematically, 3 classes of theories discussed (with complex fermions, pseu- 
doreal fermions and real fermions) are described in terms of 3 classic stochastic en- 



6 The nice result (6.7) is specific for the theories with adjoint Majorana fermions. To define in 



the Euclidean space a theory like the standard model involving chiral fermions in the fundamental 



representation is more difficult 1 1 7 
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sembles: Gaussian unitary, Gaussian orthogonal, and Gaussian simplectic. It is note- 
worthy that, though the sum rules in all 3 cases are the same for Nf = 1 [where no 
spontaneous breaking of chiral symmetry occurs, no Goldstone bosons appear, and 
the partition function has the extensive form ( |3.1| )1 , the microscopic spectral densi- 
ties and the correlators are not. That elucidates again the fact that , when deriving 
microscopic spectral densities in the stochastic matrix technique, an extra assumption 
( |5.4|) beyond using just the symmetry properties of the theory is neccessary to adopt. 
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Figure captions. 

Fig. 1. 2-point correlator, (a) Quark representation. The loop is evaluated in a 
gauge-field background to be averaged over afterwards, (b) Quasi-Goldstone contri- 
bution. 

Fig. 2. Instanton calculations for the spectral density n{X) = p(A) for different Nf. 
A is measured in units of Aqcd, and the area below the curve is normalized to 1. 

Fig. 3. Microscopic spectral density as a function of z — HVX/n. The normalization 
p(0) = 1 is chosen. Dashed lines correspond to the theoretical result Q5.7D - Full lines 
present the results of numerical calculations in the instanton liquid model. 
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